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^ , Abstract 

Grover recently presented the fixed-point search algorithm. In this letter, we 
study the fixed-point search algorithm obtained by replacing equal phase shifts 
flj I of 7r/3 by different phase shifts. 
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1 Introduction 



(N 
> 

p^ ' Grover's search algorithm and fixed-point search algorithm are both used to 

\^ . find a desired state from an unsorted database. His original algorithm consists 

f^ ' of inversion of the amplitude in the desired state and inversion-about-average 

''si" , operation[T]. In [5], Grover presented a general algorithm: Q = —I^U^^ItU, 

^^ ' where U is any unitary operation, U~^ is the adjoint of U, Ij — I—2\"/) (7I, Ir = 

I—2\t){t\, I7) is an initial state and \t) is a desired state. When U^^ = U ^ W, 
where W is the Walsh-Hadamard transformation, and I7 = |0), the general 
algorithm becomes the original algorithm. Grover showed that the desired state 
can be found with certainty after 0{VN) applications of Q to the initial state 
S . It)- Long extended Grover's algorithm[3]. Long's algorithm is Q = —^7C^~^-frC^, 

c3 ; where 1^=1- (-e*^ + 1)|7)(7|, Ir = I - {-e"^ + 1)\t){t\. When 6* = (^ = tt, 

^ . Long's algorithm becomes Grover's general algorithm. Li et al. proposed that 

**-* ' C/~^ in Long's algorithm can be replaced by any unitary operation y[4][5]. In 

[6], Galindo et al. gave a family of Grover's quantum searching algorithms. In 
0, Pang et al. generalized Grover's algorithm and applied their algorithm to 
, ; , image compression. 

C^ ' Grover's quantum search algorithm can be considered as a rotation of the 

state vectors in two-dimensional Hilbert space generated by the initial (s) and 
target (t) vectors[2]. The amplitude of the target state increases monotonically 
towards its maximum and decreases monotonically after reaching the maximum 
[S]. As mentioned in [H], unless we stop as soon as it reaches the target state, it 
will drift away. However, the number of iteration steps to find the target state 
is not an integer [5j. It means that we either stop near the target state or we 
drift away. 
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Grover presented the new algorithm by replacing the selective inversions by 
selective phase shifts of tt/S [8] . The new algorithm converges to the target state 
irrespective of the number of iterations. In [S], an algorithm for obtaining fixed 
points in iterative quantum transformations was given and it was illustrated 
that the algorithm has much better average-case behavior. In [10], Boyer et al. 
described an algorithm that succeeds with probability approaching to 1. 

The quantum search algorithm with different phase rotation angles were 
discussed in [3] [3] [5] [TT] [T^]. In this letter, we study the fixed-point search 
algorithm obtained by replacing equal phase shifts of 7r/3 by different phase 
shifts and show that the deviation for different phase shifts is smaller than for 
equal phase shifts. However, the smallest average deviation does not occur at 
different phase shifts. For the definition of the average deviation, see section 
4.1 of this letter. It is well known that the smaller the deviation is, the more 
rapidly the algorithm converges to the desired item. Let Uts be the amplitude of 
reaching the target state \t) by applying U to the start state \s) and || Uts |P= 

1 — e. For any range (/3, a) of e, we argue that a and /3 determine the phase 
shifts at which the smallest average deviation occurs. 

This paper is organized as follows. In section 2, we introduce the fixed-point 
search algorithm with different phase shifts and derive the deviation. Section 3 
is used to find different phase shifts for small deviation. Section 3 is devoted to 
study the smallest average deviation. 

2 The algorithm with different phase shifts 

2.1 Replacing equal phase shifts of 7r/3 by different phase 
shifts 

Grover presented the fixed-point search algorithm by replacing the selective 
inversions by selective phase shifts of 7r/3[8j. He described the transformation 
URsU^RtU applied to the start state \s), where 

i?. =/-(l-e"/3)|s)(s|, 
i?,=/-(l-e-/3)|t)(i|, (1) 

t stands for the target state. 

We consider the transformation URsU^RtU applied to the start state \s), 
where 



Rs = I-{l-e'n\s){s\, 

Rt^I-{l-e'')\t){t\. (2) 

6 and (p are called selective phase shifts. It is enough to consider 6 and ip to be 
in [0, tt]. When = ip = 7r/3, one recovers Grover's fixed-point search algorithm 



above. In [T3] , the fixed-point search algorithm with two equal phase shifts was 
discussed. 

We compute URsU+RtU\s) ^ (e"^ + (1 - e'-'^)(l - e^") || Uts \?)U\s) - (1 - 
e*^)t/ts|i). By the definition of the deviation in [5|, let D{9^ (p) be the deviation 
of this superposition from the state \t). Then we can derive 



Die,ip) - (1- II Uts f) II e'^ + (1 - e^m^ ~ e'*') || [/*, ff . (3) 

2.2 Simplifying the deviation for different phase shifts 

Grover let ||[/ts|| = 1 — e, where < e < 1. As indicated in ',2., \\Uts\\ is very 



small and almost l/\/N, where N is the size of the database. For example, 
e = 1 — 1/A^ > I when A^ > 4. Therefore e is close to 1. 

Substituting || Uts |P= 1 ^ e and reducing the second term of D{6, ip) in ([3|), 
we obtain 



e(l - 8(1 - e) sin(6'/2) sin(^/2) cos((6i - ip)/2) + 16(1 - e)^ sin2(6l/2) sin2(^/2)). 

(4) 



See Fig.l for D{9, ip). When 6 ^ ip, ^ becomes 

D{9, 9) = e(4(l - e) sm^{9/2) - 1)^. (5) 

When 9 = Lp = tt/3 are chosen as phase shifts, the deviation D{tt/3, tt/S) = 

em- 

2.3 The deviation does not vanish at different phase shifts 

D{9, If) in ([4]) is rewritten as follows. 

D{9, Lp) = e[(4(l - e) sin(6l/2) sin((p/2) - cos((e' - ip)/2))'^ + sin2((e' - p)l2)]. 

(6) 

From ([6|), it is straightforward that D{9,ip) = if and only ii 9 = ip and 
cos 9 = 1 — „f}_ s , where e < 3/4. It means that the deviation vanishes only 
when two phase shifts are equal. Therefore if e is known and e < 3/4, then we 
choose 9 = arccos(l — 2(i-~a ) ^^ ^^'^ equal phase shifts to make the deviation 
vanish. However, as pointed out before, e is always close to 1. 



3 Different phase shifts for small deviation 

From (HI) and ^ we calculate 

D{9,ip)-D{0,0) = 

8e(l - e) sin(6l/2) sin((6i - tp)/2)[ecos{if/2) + (1 - e) cos((26' + ip)/2)]. 

(7) 

Let us reduce (O as follows. ecos{ip/2) + (1 - e)cos((26' + f)/2) in ^= 
e(cos((y9/2) - cos((26' + ip)/2)) + cos((26' + ip)/2) = 2e sin((6l + ip)/2) sin(6l/2) + 
cos((26' + f)/2). Then ([7]) is rewritten as follows. 

D{e,ip)-D{0,0) = 

8e(l - e) sin(6l/2) sin((6' - ip)/2)[2e sin((6l + ip)/2) sin(6i/2) + cos((26l + ip)/2)]. 

(8) 

Following ([5]), we have the following results. 

Result 1. 

The deviation for different phase shifts and (p is smaller than the de- 
viation for equal phase shifts 0, i.e., D{d,Lp) < D{0,0), if < < (p and 
e > -cos((26' + ip)/2)/{2sin{{9 + ip) / 2) sin{0/ 2)) or < ip < and e < 
- cos((26' + ip)/2)/{2 sm{{9 + Lp)/2) sin(6l/2)). 

Resuh 2. 

D{e, ip) < D{0, 0) whenever < 6i < (^ and < (261 + (p) < tt. 

The following example follows result 1 immediately. 
Example 1. When = 7r/3 and 7r/3 < </? < tt, D{TT/'i,p) < D{'k/S,tt/3) for 
e > - cos(7r/3 + (^/2)/ sin(7r/6 + ip/2). 

It means that the deviation for one phase shift of n/3 and one larger phase 
shift is smaller than for two equal phase shifts of 7r/3. 

Example 2. D{0, 7r/2) < 0(9, 9) when < 61 < 7r/2 and e is large. 

The proof is as follows. 

When < 61 < 7r/4, from resuh 2 D{0,Tr/2) < D{0,0) for any e in (0,1). 
When 7r/4 < 61 < 7r/2, from result 1 D{0,tt/2) < D{9,9) for e > (sin 61 - 
cos6')/(2(sin(6l/2) + cos(6'/2))sin(6l/2)). 

This example can also be verified by computing D(9, t:/2) and reducing 
D{0,tt/2) - D{9,9) as follows. 

D{0,'k/2)-D{9,9) =4e(l-e)sin(6'/2)(sin(6'/2)-cos(6'/2)+2(l-e)sin(6l/2)cos6l). 
It is not hard to see that D{0, it/2) < D{9, 6*) if e > l-[(cos(6'/2)-sin(6'/2))/(2 sin(6'/2) cos ( 



4 The phase shifts for the smallest average de- 
viation 

4.1 The definition of the average deviation 

Usually, people pay much attention to the average-case behavior of an algorithm 
besides the worst-case behavior. For Grover's fixed-point search, Tulsi et al. 
studied the average-case behavior in [3]. As indicated in [2], \\Uts\\ is very 
small, i.e., e is very large and close to 1. Unfortunately, we don't know the 
exact value of e. However, it may be possible to know an range of e. From ((4]), 
it is clear that the deviation is a function of e. Let e lie in the range (/3, a). 
Intuitively, the deviation varies as e does in (/3, a). What is the average value of 
the deviation? It seems essential and significant to define the average deviation 
in the present letter. 

Assume that e is in the range (/3,a), where < f3 < a < 1. In terms 
of mean-value theorem for integrals, the average value D{9, Lp) over the range 
(/3, a) of e of deviation D{9, (p) is defined and calculated as follows. 



a- [3 Jp 
' -{c? - i?) + Asin(6i/2) sin((/3/2) cos((6i - (^)/2) + B sin2(6l/2) sin2(^/2)], 



a-/3'2' 



where A = -|(a2(3 - 2a) - /32(3 - 2/3)) and B = -|((1 - af{ia + 1) - (1 - 
/3)3(3/3+l)). 

We argue A < 0, S > and AjB < -1/2 in Appendix A. 

For example, when 9 — (p — n/3, the average deviation _D(7r/3,7r/3) = 
(a + /3)(a2+/32)/4. 

4.2 The phase shifts for the smaUest average deviation 

Apparently, the average deviation D{9, ip) in ([5]) is a function of phase shifts 9 
and if. It is natural to ask what phase shifts can make the average deviation 
D{9, if) the smallest. For this purpose, we need to find the minimum of the 
average deviation D{9, ip). 

To find the extremes of the average deviation D{9, cp), we compute the fol- 
lowing partial derivatives: 

dD{9,ip)/d9 ^ l/(2(a-/3))sin(^/2)(Acos((2e'-(p)/2)+Bsin6lsin(v3/2)), 

(10) 

dD{9,ip)/dip = l/{2{a- f3))sm{9/2){Acoa{{2ip~9)/2) + Bsmipsm{9/2)). 

(11) 



(9) 



The extremes of the average deviation D{6,(p) are (a + /3)/2 at the ex- 
treme points of D{e,ip): {0,0},{0,7r} and {7r,0}. Let Acos{{29 - ip)/2) + 
Bsmesm{ip/2) = in dTUD and Acos((2(p- 6')/2) + Bsin(/?sin(6'/2) = in (fTTj) . 
Then we derive the extreme points 9 — ip. Letting 9 — ip, D{9, ip) becomes 



D{9,9) = {{a^ ~ /3'^)/2 + Asm'^{e/2) + Bfim'^{9/2))/{a- P). (12) 

After rewriting, 

Let us prove that the smallest average deviation occurs at 6 = cj). 

4.2.1 The smallest average deviation for a + /3 > 1 

For large e, the range (/3, a) of e may satisfy a + /3 > 1. When a + (3 > 1, 
A/B < —1. See (3) in Appendix A. In (5) of Appendix A, we demonstrate 
when a + P > 1, the average deviation D{9, (p) < {a + (3)/2. It implies that the 
smallest average deviation occurs at equal phase shifts. From ([T^. let us find 
the smallest average deviation as follows. 

Case 1. -2 < A/B < -1 

When 9 — arccos(l + A/B) is chosen as equal phase shifts, D{9,9) reaches 
its minimum {a + P)/2 — A^ /{AB{a — /3)), which is also the minimum oi D{9, (p). 
So, arccos(l-|-A/_B), which is in [tt/2, tt], is called the smallest average deviation 
point. 

Example 3. Let /3 = and a ^ I. Then A = -4/3, B = 4/3 and 1 + A/B = 
0. D{9, 9) is calculated as D{9, 9) = cos^ 9/2> + 1/6. Straightforwardly, tt/2 is 
the smallest average deviation point at which D{9, 9) reaches its minimum 1/6. 
See Fig. 2. 

Case 2. A/B < -2 



D(9, 9) decreases as 9 increases from to tt and reaches its minimum ^ 



P 



4 (a — 1 + /3) (a^ — a + /3'^ — /3) at 6* = tt, which is also the minimum of D{ 

Example 4. Let /? = 3/4 and a = 1. Then A = -5/24, B = 13/192, 
1 + A/B = -^. D{9,9) becomes (13/192)(27/13 + cos6l)2 + 73/312 and Z)(6l, 6*) 
reaches its minimum ^ at = tt. See Fig. 3. 

4.2.2 The smallest average deviation for a + j3 <\ 

When a + /3 < 1, —1 < A/B < —1/2. Then it can be shown that when 
9 = arccos(l + A/B), D{9, 9) reaches its minimum (a + (3)/2 - A^/{'iB{a - (3)), 
which is less than (a + /3)/2, and when 9 = n, D(9,9) reaches its maximum 
^ + A{a - 1 + l3) {a'^ - a + f3'^ - 13), which is greater than (a + /3)/2. 

As stated already, the extremes of D{9,(t)) at {0, 0},{0, tt} and {7r,0} are 
{a + P)/2. Consequently, when a + f3 < 1, the smallest average deviation occurs 



at 6* = = arccos(l + A/B), which is in (tt/S, 7r/2), and the maximal average 
deviation occurs at 6 = (j) = n. 

Remark: 

Conclusively, for any range (/3, a) of e, a and (3 determine the phase shifts 
at which the smallest average deviation occurs, the smallest average deviation 
points are greater than 7r/3 and the smallest average deviations are smaller than 
the average deviation Z)(7r/3,7r/3) for equal phase shifts of 7r/3. 



5 Summary 

In this letter, we demonstrate the possibility of the fixed-point quantum search 
algorithm with two different phase shifts. Intuitively, not only there are more 
choices for phase shifts to adjust an algorithm for future physical realization, 
but also we can find some different phase shifts for small deviation. Thus, more 
loose constraint opens a door for more feasible or robust realization. In this 
letter. We also show that the smallest average deviation can be obtained by 
choosing the following equal phase shifts. Let (/3, a) be the range of e. Then if 
A/B > —2, then arccos(l + A/B) is chosen as equal phase shifts. Otherwise, 
the closer to tt the equal phase shifts are, the better. 
Appendix A 

1. Theproof of B > 

We can factor B as follows. 

B = -|((1 - a)3(3a + 1) - (1 - /3)3(3/3 + 1)) = | (a - /3) D, where D = 
3a^ - 8a^ "+ Sfia^ + 6a - 8/3a + Sfi'^a + 6/3 - 8/3^ + 3/3^. So, we need only to 
show that D > 0. D can be rewritten as (a + /3)[3(a + (3)^ - 8{a + /3) + 6] - 
a/3[6(a + /3)-8]. There are two cases. Case 1: 6(a + /3)-8<0. Clearly D > 
since 3(a + /3)2 _ 8(a + /3) + 6 > 0. Case 2: 6(a + /3) - 8 > 0. Also, D>0 since 
(a + /3)[3(a + /3)2-8(a + /3) + 6] > (a + /3)2[6(a + /3) - 8]/4 > a/3[6(a + /3) - 8]. 

The proof was given by Mr. P.Y. Sun. 

2. The proof of A/B < -1/2 

We demonstrate i? > in Appendix A. Hence, it is easy to see A/B < —1/2 
if and only if 2A + B < 0. By factoring, 2A + B = ^ {a - P) E, where E = 
3a^ - ^o? + 3/3a2 - 4/3a + 3/3^0: - 4/3^ + 3/3^. We only need to argue -B < 0. 
Letting a = /3 + 7, where < 7 < 1, 

E = 12/33 _^ Y%0^-i + 12/372 + 'i-f - 12/32 _ ;^2/37 - 472 

= (12/?3 + 12/3^7 - 12/32) j^ (g^2^ ^ ^^2/372 - 12/37) + (37^ - 472) = 

= 12/32 (/3 + 7 - 1) + 6/37(/3 + 27 - 2) + 72(37 - 4) 

= 12/32(a - 1) + 6/37(a + 7 - 2) + 72(37 - 4) < 0. 

3. The proof of A/B < -1 when {a + (3) > 1 

Since A + B = 4{a - 1 + l3) {a - l3) {a"^ - a + 13'^ - l3), clearly A/B < -1 
when {a + (3) > 1. 

4. The proof of yl < 



By factoring, we obtain A = ^{a ~ f3)C, where C — 2(a^ + a/3 + /3^) — 
3(a + f3). By letting a = /3 + 7, where < 7 < 1, C = 6a/3 + 27^ - (6/3 + 87). 
It is easy to know C < since 6a(3 < 6/3 and 27^ < 87. 

5. The proof of D{e, if) < {a + /3)/2 when a + /3 > 1 

When a + /3 > 1, we show that A/B < — 1 in (3) of this appendix. To show 
D{0, f) < (a + /3)/2, from © we only need to show Asin(6'/2) sin((^/2) cos((6'- 
(p)/2) + Bsm'^{6/2) sin^((p/2) < 0. Let us argue this as follows. 

Asin(6»/2) sm{(p/2) cos((6l - ip)/2) + B sin^{e/2.) sin^{ifi/2) 

= i?sin(0/2)sin(^/2)cos((0 ^ ^)/2)[4 + Si^WW^j. 

Mpvt ^p rnmniitp ^ + sin(e/2) sin(y/2) _ ,4 , 1 1 cos((e+y)/2) c:„„p A <- 

iNext we compute -g + cos((e-y)/2) " s + 2 ~ 2 cos((e-i^)/2) • ^™ce 5 S 

_-. A I 1 _ 1 cos((9+y)/2) ^ _1 _ 1 cos((6)+y)/2) _ _!/-, cos((g+y)/2) x _ 

-•-'5^2 2 cos((6(-i^)/2) — 2 2 cos((e-¥3)/2) 2^-^ ^ cos((e->^)/2) '' 

cos(e/2) cos((/)/2) ^ p. 
cos((9-v)/2) ^tJ- 
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